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Abstract 

The general theory on exact boundary controUabihty for general first order quasilin- 
ear hyperbolic systems requires that the characteristic speeds of system do not vanish. 
This paper deals with exact boundary controllability, when this is not the case. Some 
important models are also shown as applications of the main result. The strategy uses 
the return method, which allows in certain situations to recover non zero characteristic 
speeds. 
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1 Introduction and main results 

The general theory on exact boundary controllability for general first order quasilinear 



hyperbolic systems requires that the system has non vanishing characteristic speeds |19, 2C]. 
Several papers have dealt with hyperbolic systems having a vanishing or an identically zero 
characteristic speed, under various assumptions. For systems with identically zero charac- 
teristic speeds, a general result on exact controllability has been obtained by using internal 
controls |21]. It is also possible to get in this case partial controllability by boundary controls. 



if some eigenvalue of the system is equal to zero identically [^] . A steady state controllability 
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holds for some special hyperbolic models with vanishing characteristic speed as Saint- Venant 



equations (or shallow water equations), see Gugat |12]. For what concerns the system of 
isentropic gas dynamics (which contains the Saint- Venant model), a more general boundary 
controllability result for (non constant) BV solutions was obtained by the second author in 

In this paper, we will discuss exact boundary controllability for a general hyperbolic 
system which admits a vanishing characteristic speed. 

Consider the following first order quasilinear hyperbolic system 
nil nil 

— + A{u)— = 0, {t,x)e[0,T]x[0,L], (1.1) 

where u = (ui, • • • , UnY^{t, x) is the state of the system in some nonempty open set il. C M" 
and the n x n matrix A belongs to (7^(0; M"^"). 

Let u* € il. he fixed. Assume that A{u*) has n real distinct eigenvalues: 

Ai(n*) < • • • < Xm-i{u*) < A„(n*) = < A^+i(n*) < • • • < A„(n*), (1.2) 

for some m G {1, • • • ,n}, which are the characteristic speeds at which the system propagates. 
Thus in a neighborhood of the equilibrium u = u*, the system is strictly hyperbolic and A{u) 
has a complete set of left (resp. right) eigenvectors li{u), ■ ■ ■ ,ln{u) (resp. ri{u), ■ ■ ■ ,r„(u)): 

li{u)A{u) = \i{u)li{u) {resp. A{u)ri{u) = Xi{u)ri{u)), i = l,---,n. (1.3) 

Without loss of generality, let us assume that 

li{u)rj{u) = 6ij, i,j = l,--- ,n, (1.4) 

where 6ij is Kronecker's symbol. Reducing 17 if necessary, we assume that 

Vj G {1, . . . ,771 - 1}, Xj{u) < and Vj G {m -M, . . . ,n}, Xj{u) > 0, Wu £ Q. (1.5) 

Now the question is: is it possible to realize the local exact controllability near the equilibrium 
u = u* only by using boundary controls? 

In order to overcome the difficulty of a characteristic speed vanishing at u* , we assume 
the following hypothesis: 

(H): for all e > 0, there exists a = (ai, • • • , am,_i, Om+i, • • • , On) G L°°{0, 1; R"~^) with 

II«IIl°°{0,1;R"-1) < (1-6) 

such that the solution z £ C'^([0, 1];M'") of the ordinary differential equation 

-=^a,{s)r,{z), z{0)=u*, (1.7) 

satisfies 

A„,(z(l)) / 0. (1.8) 
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The main result of this paper is the following theorem: 

Theorem 1.1. Let (1.2) and (H) be true. Then, for any 6 > 0, there exist T > and u > 



such that, for all ^p,tp £ C^([0, L];]R") satisfying 

IIV(-) -""llcHIO.L]) < IIV'(-) - ^*llci([0,L]) < ^' (1-9) 

there exists u G C^([0,r] x [0,L];M") such that 

(111 Cfll 

— + A{u)— = 0, V(t,x)G[0,r]x[0,L], (1.10) 

u{0,x) = ip{x), VxG[0,L], (1.11) 
u{T,x) =i{;{x), VxG[0,L], (1.12) 
h(ir)-n*|bi([o,L]) <'5, VtG[0,r]. (1.13) 

The hypothesis (H) seems quite difficult to check. However, we have some sufficient 
conditions of (H) relying on Lie brackets. 

Proposition 1.1. The following properties are sufficient conditions for (H) to hold: 
(HI): there exists j G {1, • • • ,n}\ {m} such that VXmiu*) ■ fj{u*) ^ 0, 
(H2): there exist j, k £ {1, ■■ ■ ,n}\{m} such that VXm{u*) ■ [rj,rk]{u*) ^ 0, 
(H3): A G C°°(J7;M"^") and there exists h G Lie{ri,--- ,rm-i,rm+i, ■ " ^^n}, such that 

(H4): A G (:;°°(17;M"^") and {/i(n*), h G Lie{ri,--- , r^_i, r„+i, • • • ,r„}} = and u* is 
in the closure of {u £ Q : \m{u) 7^ 0}. 

Here Lie{ri,--- , r^-i, Tm+i, • ■ ■ ifn] denotes the Lie algebra generated by the smooth 
vector fields n,. . . ,rm_i,rm+i,. . . ,rn. 



Proof of Proposition Ll . It is a consequence of Chow and Rashevski's connectivity Theorem 
(see for instance Theorem 3.19, p. 135]) that (H4) implies (H). Next we notice that both 
(HI) and (H2) clearly imply (H3). So we have left to prove that (H3) implies (H). From (H3) 
we deduce that the exists a direction b G M" obtained by p successive Lie brackets and such 
that VXm{u*) • /i / 0. We use [15, Lemma 1, p. 456] to deduce that there are controls a 



which are arbitrarily small in L°° norm such that the corresponding solution of (|1.7] ) satisfies 
^(4Pti/(i+p)) = u* + tb + o{t) as t ^ 0. The conclusion follows. □ 



Remark 1.1. Theorem 1.1 can be regarded as a local boundary controllability result because 
one can drive any initial data ip to any desired data ip near u = u* without using any internal 
controls. However, since the characteristic speed Xm may change its sign during the control 
period, it is difficult to describe the exact distribution of boundary controls. To overcome this 
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difficulty, we consider the system without boundary conditions (which is consequently under- 
determined), and aim at finding the solution u itself. In the conservative case (where A(u) 
is a Jacobian matrix Df{u)), the solution that we determine can enter the general theory of 
initial-boundary problems for systems of conservation laws, see in particular Amadori Q and 
Amadori and Colombo Q. 

Remark 1.2. Due to the hypothesis (H), one can drive the possible vanishing characteristic 
speed Am to be nonzero after sufficiently long time by only using boundary controls. However, 
if some characteristic speeds of the system are identically zero, the approach of this paper is 
not valid anymore. Is boundary controllability possible in such cases, even for some special 
models? Up to our knowledge, this question remains open. 

Remark 1.3. We could treat the case where A £ (7^(0; M"^"), see in particular Remark 
P) below. 



The main idea to prove Theorem 1.1 is to use a constructive approach and the return 
method In our framework the method consists in constructing a trajectory n G C^([0, T] x 
[0, L]; M") of the system ([LI]), close to u* such that 

u{0,x) =u{T,x) = u*, VxG[0,L], (1.14) 

and that the linearized equation around u is controllable. Note indeed that the linearized 
equation around u* is not controllable. Based on this, we can construct a solution u G 
C^([0,T] X [0, L];]R") to the system ( |1.1| ) which connects the initial and final data (which 
have to be sufficiently close to u*). 

As a matter of fact, we will quite not use the linearized equation. Instead, we use an 
argument of perturbation of the trajectory u and then reduce the original control problem to 
a boundary control problem without vanishing characteristic speeds, which has been solved 
by Li and Rao [^]. In the framework of systems of conservation laws, the return method has 
also been used in |8|, [l^, 17 1, see also |^]. For other applications of the return method, see 
1^] and the references therein. 

Without loss of generality, we may assume the equilibrium u* to be 0, replacing u by 



u — u* as the unknown in the system (LI) if necessary. For the convenience of statement, we 
denote by C various positive constants in the whole paper which may change from one line 
to another. 

The organization of this paper is as follows: in Section 2 we construct the special trajectory 
u G C^([0, T] X [0, L]; E") of the system (lO) which starts at and returns to 0, and such that 



the equation linearized around u is controllable. Then we prove the main result, Theorem 1.1 



in Section 3. Some important applications are shown in Section 4, including Saint- Venant 
equations (shallow water equations), 1-D isentropic gas dynamics equations, 1-D full gas 
dynamics equations and Aw-Rascle model on traffic flow and its generalization. Finally in 
Appendix A, we establish a technical result. 
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2 Construction of the trajectory u 

Definition 2.1. Let j G {1, • • • ,n} and G Q. Let s G [— eo,eo] ^ Uj{s) G O be the orbit 
of the eigenvector field rj starting at (or rarefaction curves): 

^ = r,(C/,), [7,(0) = nO, (2.1) 

where Eq > is a small constant. Let ^j{s, •) be the corresponding flow map when s varies, 
i.e., 

$,(s,nO) :=[/,(s), VsG[eo,eo]. (2.2) 

Remark 2.1. For all s G [— eoi^o] one has n"*" = ^j{s,u~) <^=^ u~ = ^j{—s,u~^). 

Our first proposition concerns simple waves which one can use to modify the state in 
[0,L]. 

Proposition 2.1. Let j G {1, • • • ,n} \ {m} and 

There exist C > and Eq > 0, such that for all e G (0, £q\, all u~ ,u~^ £ il. satisfying 

I"""!, I^^"*"! < e and = $j(s, u^) for some s such that \s\ < e, (2-4) 
there exists u G C^{[0,T] x ]R;]R") such that 

— + A{u)— = Q, V(t,x) G [0,r] xM, (2.5) 

u{Q,x) = u~, Vx G [0,L], (2.6) 
ii(r,x) = n+, VxG[0,L], (2.7) 
\\u{t,-)\\cHR)<Ce, VtG[0,T]. (2.8) 

Proof: Without loss of generality, we may assume that j G {1, • • • ,m — 1} (the case where 
j G {m + 1, • • • , n} can be treated similarly by symmetry in x, that is, replacing x hy L — x 
if necessary). 

In view of (|1.2D and (|2.3| ), there exist ei > and 77 > small enough such that 

T > max (2.9) 

\u\<ei \Xj{u)\ 



Let e G (0,ei] and n G be such that ( |2.4| ) holds. By Definition 2.1, it is easy to 
see that 

\^j{s,u-)\ <Ce, Vs G [-|s|,|s|]. (2.10) 
Let P G C^((0, 1);M) be such that 

^(0)^61 = 1. (2.11) 



Then we let 



/3(e) :=-/?(-), (2.12) 
7] rj 



which gives that /3 G C^f ((0, r/); M) and 



V 



f3{e)de = s. (2.13) 







From the above, the ordinary differential equation 

^ = mrM^ y(0) = n-, (2.14) 
admits a unique solution y(-) = <^j{a{-),u~) G C2([0, ??]; M"), where 

cj(s) := /" p{e)de, Vse[0,ry]. (2.15) 

Let 

III", X < L, 

y{x-L), L<x<L + r], (2.16) 
U"*", X > L + 7]. 

In the following, we will denote by C'^ (M) the space of functions of class whose derivatives 
up to order k are bounded on M (and the norm || • ||cfc(K) is in fact the 
Then by (|t|), (|2T2|) , (|2T|) and (|1^), we obtain that 



\f\\co(R) ■= snp\ip{x)\ < Ce, (2.17) 
|¥.'||cO(M) := sup \ip'{x)\ <C-< Ce. (2.18) 



Now we focus on the Cauchy problem of ( |2.5D on R with the initial condition 

u{0,x) = ip{x), VxGM. (2.19) 
It is classical that there exists a unique solution to the Cauchy problem ( |2.5| ) and ( p. 19] ) 



in small time; see for instance |16, p. 55]. Let us prove that: for the fixed time T > 0, 
if e is sufficiently small, the Cauchy problem ( |2.5[ ), ( |2.19| ) admits a unique solution u G 
C2([0,r] X M;M") such that (^ to (|2^ hold. 

To show that, it suffices to obtain a uniform a priori estimate of the solution in (see 
|16, Theorem 4.2.5, p. 55]). In order to obtain such an a priori estimate, we assume that 



the Cauchy problem (p. 51) , (|2.19|) admits already a solution u G C^([0,ro] x M;M") for some 

To G (o,r). 

For any i G {1, • • • ,n} and any point {t,x) G [0,To] x M, we can define the i — th 
characteristic curve ^ = ^j(r) passing through {t,x) by 

(if 

^ = X,{u{t,0), m = ^- (2.20) 
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Introducing 



I.e., 



Vi := k{u)u, Wi := li{u)^, i = 1, • • • , n, 



we know that Vi,Wi {i = 1, ■ ■ ■ ,n) satisfy the following (see [16, p. 47ff] and [[LSl]): 



-T^ = ^Pikiiu)vkWi, i = l,- 



dit 

dw. 



-77 = ^liki{u)wkWi, i = 1, 

dit 



• • • ,n, 



k,l 



where 



d _ d X ( \ ^ 



(2.21) 
(2.22) 

(2.23) 
(2.24) 

(2.25) 



denotes the derivative along the i-th characteristic, and where Pikh^iki S C^{i^;W^) satisfy 
in particular 



'yikk{u) = 0, V«,fc e {!,••• A; / i, 
-yuiiu) = -VXi{u)ri{u), i = l,---,n. 



By (|2.26D - (^l27|) , (|2.24D can be written as 

dw, 



dd 



'^Jiki{u)wkWi - {V Xi{u)ri{u))wf , i = l,' 



, n. 



(2.26) 
(2.27) 



(2.28) 



ki^l 



Combining (p^)-( [2l8|) and (2.22), noticing (0), we have 



\vi{Q,x)\ = \li{Lp{x))Lp{x)\ < Ce, \/x e R,\/i G {!,••• ,n}, 

Wi{o,x) = k{ip{x))ip'{x) = 0, yx eR,yie {!,■■■ ,n}\{j}, 

\wj{0,x)\ = \lj{ip{x))ip'{x)\ < Ce, Vx G M. 
As in the proof of [|l^, Theorem 4.2.5, p. 55], we first assume that 

\vi{t,x)\<l, \wi{t,x)\<l, V(t,x) G [0, To] X M,Vi G {!,••• ,n}. 
By (2.28) when i ^ j and ( p. 30 ) , we deduce that 

w^{t,x)=0, y{t,x) G [0,To] xM, Vi G {!,••• ,n}\{j}, 



which then reduces ( |2.28| ) when i = j to 

dwi 



it 



(2.29) 
(2.30) 
(2.31) 



(2.32) 



(2.33) 



(2.34) 
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By comparing the norm of the solution ( |2.23| ) and ( 2.34| ) to the solution of the ordinary 
differential equation x = ||VAj.rj||x^, noticing also ( 2.29| ) and ( 2.31| ), we deduce that the 
following estimates hold 

\vi{t,x)\ <Ce, V(i,x) G [0, To] X M, Vi G {!,••• (2.35) 
\wj{t,x)\ < Ce, y{t,x) G [0,ro] X M. (2.36) 

Combining ( 2.33| ) and ( 2.35[ )-( ^l3^ ), there exists Eq G (0,ei] small enough such that for any 
e G (0, eo] the assumption ( |2.32| ) is indeed satisfied, and the uniform a priori estimate 

\\u{t,-)\\cHR)<Ce, VtG[0,To] (2.37) 

holds for all Tq G (0, T). This proves the existence of the solution u G C^([0, T] x M; R") (see 



again |16, Theorem 4.2.5, p. 55]). Moreover, since ( p.9D implies 

T > ^7^, (2.38) 
\Xj{u+)\ 

we derive from the fact 

Wi{0,x)=0, Vx G [L + r/,c3o),Vi G {!,••• ,n} (2.39) 

that 

Wi{T,x) = 0, Vx G [0,oo),Vi G {I,--- ,n}, (2.40) 
which in turn implies that 

u{T,x) = const. = lim u{T,x) = u~^, Vx G [0,oo). (2-41) 



This concludes the proof of Proposition 2.1. □ 



Remark 2.2. In the proof above, the information travels from right to left through the 
boundary x = L. For i G {1, . . . ,m — 1}, the information would travel from left to right 
through the boundary x = 0. 



on 



Remark 2.3. Because ( |2.5| ) is an autonomous system, the conclusion of Proposition 2.1 
[0, T] can be achieved on [to, to + T] for any to £ by translation in time. 

The next proposition prove that one can approximate the trajectory given by (|l]^) by a 
trajectory composed of simple waves. 

Proposition 2.2. There exist C > and Eq > such that the following holds. For any 
e G (0,eo], any a = (ai, • • • ,0^-1,"™+!, •••,««)£ i^°°(0, 1;R"~^) satisfying 

II"IIl°°(o,1;R"-i) < (2.42) 

we consider z G ^^([0, 1];M") the solution to the ordinary differential equation 

dz 

- = Y^a,{s)r,{z), z(0) = 0, (2.43) 



Then, for any r] > 0, there exist p G N, ii, • • • ,ip £ {1, • • • ,n}\ {m} and ti, • • • , tp S M such 
that 



^\ti\<Ce, (2.44) 



1=1 

|z(l) - •) o • • • o $,,(t2, 0) (ti, 0))| < 7]. (2.45) 

Proposition will be established in Appendix A. The next proposition, which establishes 
the existence of the special trajectory u, is the principal of this section. 

Proposition 2.3. Let K be a compact subset of ft. There exist C > and > such that 
the following holds. For any e £ (0,eo], there exist T > and a state u* £ K satisfying 

Ai(¥*)/0, Vi G {!,••• ,n}, (2.46) 

there exist p £ N and times = tq < ri < • • • < r2p+i = T with 

Tp+i -Tp> max f , (2.47) 

1=1. ..n \Xi[U )\ 

and a function u £ L°°{{0,T) x ]R;IR") such that 

ui[r,_„n]xR^C\[Ti_,,Ti]xR;R^), V/G{l,...,2p+1}, (2.48) 

^|[o,T]x[o,L] e C\[0,T] X [0,L];E"), (2.49) 

du cPu, 

— + Aiu)— = for it, x) in each [n-i^n] x M, Ml £ {1, . . . , 2p + 1}, (2.50) 
ot ox 

u{0,x) =u{T,x) =0, VxG[0,L], 

u{t,x) = U*, \/t £ [Tp,Tp+i],\/x £ [0,L], 

\\u{t,-)\\c^R) < Ce, fort in each [ti_i,ti], \/l £ {l,...,2p + l}. (2.51) 

Proof: By Proposition and the hypothesis (H), we can deduce that there exist C > 
and ei > such that for any e £ (0, ei], one can find p € N and ii, - ■ ■ ,ip £ {1, • • • ,n'} \ {m}, 
ti , • • • ,tp £M. such that ( |2.44 ) applies and 



Xm{<^^,{tp, > • • • o ^,^{t2, •) O ($i^(ti,0))) / 0. (2.52) 

And thus 

X^{^^,{tp,■)o■■■o^^2{t2,■)°{'^^^{tl,0))) ^0, Vi G {1, • • • , n}. (2.53) 

We let 

U* := ^i^{tp, •) o • • • o $,^(t2, •) o ($,^(ti, 0)). 
Now for every / £ {1, • • • ,p}, let 
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and in addition 



I 

n := ^Tfc for / = 1, . . . ,p, 

k=l 

L 

Tp+i := Tp + max + 1 

^ ^ j=l...n |Aj(M*)| 

p 



(2.55) 



T; := Tp+i + ^ Tfc for / = p + 2, . . . ,2p+ 1. 

k=2p+2-l 

Observe that there is a symmetry with respect to the central time interval [rp,rp^_i], that is, 
[Tp^i,Tp] is symmetric of [Tp+i,Tp+2], etc. 



Applying Proposition 2.1 and Remark 2.2 with 

u_ = $i;_^(ti_i,-)o---o^>i2(*2,-)°(^n(il,0)) and u+ = ^>i, (t;, •) o • • • o ^.^^fe, •) ° ('J'n (^l, 0)), 

for I = 1, . . . ,p, we deduce that provided that Eq is small enough, for any e G (0,eo], there 
exists n' G C2([tj_i, r;] x M;]R") such that 



— + A(u')^ = o, y{t,x) e[Tl_i,Tl]x 

ji{Ti_i,x) = n_, Vx G [0, L], 
Ii'(r;, x) = u+, VxG[0,L], 
ll^*'(i,-)||cim) < Ce, Vt G h_i,ri]. 



T := r2p+i. 



Then, we let 
Finally, letting 

u\t,x), {t,x) G [ti-i,ti] xR,l = !,■■■ ,p, 

u{t, x) := lu*, [t, x) G [Tp, Tp+i] X M, 

u^P+^-^{Ti-t,L-x), {t,x) G h_i,r/] x]R,/ = p + 2,--- ,2p + l. 
we can see that u G L°°((0,T) x ]R;R"') satisfies the required properties. 



(2.56) 

(2.57) 
(2.58) 
(2.59) 

(2.60) 



(2.61) 



□ 



3 Proof of Theorem 

In order to conclude the proof, we will use a perturbation argument together with a result 
by Li and Rao |2^]. First, we have the following perturbation result. 

Proposition 3.1. Consider K d ^ a nonempty compact subset. Let T > 0. For any 

u G C2([0,r] X R; JO satisfying 



— + A{u)— = Q, V(t, x) G [0, r] X M, 
u{Q,x) ='4j{x) VxGM, 



(3.1) 
(3.2) 



10 



there exist uq > and C > such that for any v G (0, vq) and any ip G C^(M; Q) satisfying 

-uiO,-)\\ci(m <i^, (3.3) 



then the unique maximal solution u £ C^{[0,Tq] x M; Q) of 

nil nil 

— + A{u)— = 0, V(t,x)G[0,r]xM, (3.4) 

u{0,x) =i{;{x), VxGM, (3.5) 

is defined on [0, T] x M and satisfies 

\\u{t,-)-u{t,-)\\cHR)<Cu, ytG[0,T]. (3.6) 

Proof: Given t/j G C^(M;r2), there exists a local in time solution u G C^([0, Tq] x M) of 
j^)- (|3.5| ). We show in the same time that u does not blow up before T and that (|3.6| ) holds. 
For that, let us make the difference of (^J|) and (|3.4[), we get 



-{u-u)+A{u) — {u-u) = iA(u)-A{u)) — , V(t,rE) G [0,T] xM, (3.7) 
u{0,x) -u{0,x) ='il:{x) -i){x), VxGM. (3.8) 

By Gronwall's inequality we deduce that 

\\u{tr)-u{t,-)\\coim)<C\\^l;-ip\\coiu)<Cu. VtG[0,r]. (3.9) 

Differentiating with respect to x and observing that u is of class C^, we can use the 

same Gronwall argument to infer (^]^) and that the maximal solution is defined on [0, T]. 

□ 

Remark 3.1. We could use only a regularity assumption on u provided that this u has 



the particular structure given by Proposition 2.1. While the estimate (|3.9D should be replaced 



by a weaken one (but sufficient for the proof of Theorem ^]^) : 

\\u{T,.)-u{T,.)\\ci(^[o,L])<Cu. (3.10) 



Remark 3.2. As previously, the conclusion of Lemma on [0,T] can be achieved on 
[to, to + T] for any to £ by translation in time. 



Proof of Theorem 1.1: Again, we may assume the equilibrium u* to be 0, otherwise we 



can replace uhy u — u* as the unknown in the system (|l.lD . 



By Proposition 2.3, we can deduce that: there exist C > 0, eo > and T > 0, such that 
for any e G (0,eo], there exists u G L°°{{0,T) x R;W) such that (^l6|)-(^) hold. 
For every / G {1, • • • let t; be given by ( 2.55| ). Let 



u°{0,x) = ip{x), Vx G [0,L]. (3.11) 
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The proof relies on a induction argument on By Proposition |3.1| , we see that there exist 
C > 0, > and > 0, for any e G (0, ej and any z/ G (0, f,], if 

\\u'-\ti.i, •) - u(rz_i, Ollciao.L]) < (3.12) 
then there exists u' G C^([r/„i,ri] x R;R"') such that 

— + = 0, V(t,x) G h_i,r;] xM, (3.13) 

u\ti.i,x) = u^~\ti^i,x), Vx G [0, L], (3.14) 

IK'(i,-)||cMK) <Ge + '^'^> Vt G [Ti_i,ri], (3.15) 

||n'(T«,-)-n(rz,-)||cM[o,L]) <'^'^- (3-16) 

Therefore, there exist C > 0, ef > and > 0, such that for any e G (0, ef ] and for any 
u G (0,1/f], if 

M\c^{[QM) ^ ^' (3-17) 
then there exists G C^{[0,Tp] x [0,L];M") such that 



+ A{u^)— = 0, y{t,x) G[0,Tp]x[0,L], (3.18) 



'u*'(0,x) = (^(x), VxG[0,L], (3.19) 
h'(ir)lbi([o,L]) <Ce + Ci/, VtG[0,rp], (3.20) 
||n''(rp,-) -u{Tp,-)\\ci([o,L]) < Cu. (3.21) 



In the same way and in view of Remark 3.2, there exist C > 0, eb > and u^, > 0, such that 
for any e G (0, eb] and for any G (0, z/b], if 

llV'llcM[o,L]) < (3-22) 
then there exists G C^{[Tp+i,T] x [0,L];]R") such that 

— + A{v!")— = 0, V(t, x) G [Tp+i, T] X [0, L], (3.23) 

u^(T,x) = V(x), Vx G [0,L], (3.24) 

h''(i,-)lbi([o,L]) < C^e + C'^, VtG[rp+i,r], (3.25) 

||n''(rp+i, •) - u{Tp+i, ■)\\cH[o,L]) < Cv- (3.26) 

Now we can apply the result of Li and Rao [^] near the equilibrium of u{Tp, •) = 
u(Tp+i,-) = u* G due to ( p. 47 ) there exists Um > 0, such that for any v G (0, z^m], 
if \\u{Tp,-) — w*||(7i([o,L]) and ||u(rp+i,-) — w*||(71([o,l]) are small enough, there exists G 
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C^{[Tp,Tp+i] X [0,L];M") such that 



0, y{t,x) £[Tp,Tp+i]x[0,L], 



u'^{Tp,x) = u^{Tp,x), Va;G[0,L], 
n'"(rp+i,x) = n''(rp+i,x), Vx G [0, L] 



(3.27) 

(3.28) 
(3.29) 
(3.30) 



Combining ah of the above, there exists C > such that for any 6 > 0, there exist e > 
and > smaU enough, such that for any ip,ip £ C^([0, L];]R") satisfying 



M\CH[0,L]) ^ ^' ll^llci([0,L]) < Z^, 

one can construct u £ C^{[0,T] x [0,L];M'") by 

\Ht,x), V(t,x) G [0,Tp] X [0,L], 
u{t,x) = I u'^{t,x), V(t,x) G [Tp,rp+i] X [0,L], 
^n*'(i,x), V(t,x) G [Tp+i,T] X [0,L]. 

Now this function u clearly satisfies 

^ + A{u)^ = 0, V(t,x)G[0,T]x[0,L], 

u(0, x) = V9(x), Vx G [0, L], 
n(T, x) = ip{x), Vx G [0, L], 

Mt,-)\\cmQ,L])^C£ + Cu<6, vt G [o,r]. 



(3.31) 



(3.32) 



(3.33) 

(3.34) 
(3.35) 
(3.36) 



This finishes the proof of Theorem |L1 



4 Some models 



Model 1: Saint- Venant equations (shallow water equations) 12, 



where (7 > is the gravity constant. Let U = {H, y)*^, ( |4.4| ) is reduced to 

Ut + A{U)U, = 



with 



A{U) 



V H 

9 V 



(4.1) 



(4.2) 



(4.3) 
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By the study of Model 2 (see below), Theorem LI can be applied to (O) near the equilibrium 



U* := {H*, V*) where V* = \/gH* with H* > oi near the equilibrium U* := {H*, V*') where 
V* = -y/gH^ with H* > 0. 



Model 2: 1-D isentropic gas dynamics equations in Eulerian coordinates |11 



dp dm ^ 

dt dx ' 
dm d ,m'^ 

dt dx p 



(4.4) 



where 

p = Kp^ {K>0, 1< 7 < 3). (4.5) 

We can see (0) is a special case of ( ^ ) when p = gp^ /2. 
Moreover, let U = (p, u)*^, ( ^ ) is reduced to 

Ut + A{U)U^ = ^ (4.6) 

with 

MU) = ( ^ [ ) . (4.V) 

The characteristic speeds and the corresponding eigenvectors are 



\i{U)=u-^Mp), X2{U)=u + yMp), (4.8) 
^l([/) = (^4=,-l)*^ r2{U) = {^£=,ir- (4.9) 



Let U* := {p*,u*) where u* = \/p'{p*) with p* > 0, that is, the fluid reaches the sound 
speed. Then it is easy to check that 



Ai(C/*) = 0<A2(C/*) = 2vy(p*) (4.10) 
and the hypothesis (HI) is satisfied as: 

VAi(C/*) •r2(^7*) = > 0. (4.11) 



Similarly, if we let U* := {p*,u*) where u* = \/p'{p*) with p* > (which is the symmetric 
case of the latter) , one can see that 



Ai(C/'^) = -2vy(7)<A2(t/') = (4.12) 
and the hypothesis (HI) is satisfied as: 

VA2(?7*) • ri{U*) = < 0. (4.13) 



Therefore, Theorem 1.1 can be applied to ([4.4[) near the equilibrium U* or U*. 
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Model 3: 1-D full gas dynamics equations in Eulerian coordinates |22]: 

dp d , . 



d_ 

d \ ,v? 



(4.14) 



dt 



d_ 

dx 



,u 



pn(— + e) +pu 



0. 



Assume the gas is polytropic, so that 



and 



e = CyT = icy > 0, R>Q) 

P 



p = ke^v (A; > 0, 1 < 7 < 3). 



(4.15) 
(4.16) 



Thus, on the domain of p > 0, we have Pp > 0, Ppp > and p^ > 0. Model 3 generalizes 
Model 2 if we let m := pu and S = Sq & R. 



Let U = {p,u,S) , then (4.14) can be rewritten as 



Ut + A{U)U, = 0, 



with 





u 


P 









Pp 


u 


Pi 






p 




p 




v 








u 


) 



tr 



A{U) 



The characteristic speeds and the corresponding eigenvectors are 

Ai(;7) = n-c, \2{U) = u, \^{U) = u + c, 
rl(C/) = (p,-c,0)*^ r2{U) = {p,,G,-ppt, r^{U) = {p,c,{)) 

with c = -y/p^- 

Let U* := (/9*, 0, S*) where p* > 0, 5* G M, then it is easy to check that 

Ai([/*)<A2(t/*)=0<A3(f/*) 
and the hypothesis (HI) is satisfied as: 

VA2(f/*) -^([7*) = -c(?7*) < or VA2([/*) • r3([7*) = c(C/*) > 0. 



(4.17) 



(4.18) 



(4.19) 
(4.20) 



(4.21) 



(4.22) 



Therefore, we can apply Theorem LI to obtain boundary controllability for ( 4.14 ) near the 
equilibrium U* . 

Model 4: AR and MAR traffic flow system @, |]: 



— {p{u + pip))) + —{pu{u+p{p))) = 0, 



(4.23) 
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with 



and 



p = p"i (7>0), 



P 



'i _ 1 



(7 > 0, po > 0). 



We deduce system (AR) form system ( MAR ) by letting po = +00. 
Let U = {p,uy, dOsD is reduced to 



Ut + A{U)U, = 0, 



with 



AiU) 



(AR) 
(MAR) 



^0 u — pp'{p) J 
The characteristic speeds and the corresponding eigenvectors are 

Xi{U) = u-pp'{p), X2{U)=u, 
rl(C/) = (l,-/(p))*^ r2(C/) = (l,0)*'-. 

Let U* := {p* ,u*) where u* = p*p'{p*) with < < po> then we have 

Ai(C/*) = 0<A2(C/*)=/>V(p*)>0, 

and the hypothesis (HI) is satisfied as: 

VAi(C/*) • r2(f/*) = - P*P"{P*) = -l{p*r\-, - -)-^~\^ + 7) < 0. 

P Po Po 

Similarly, if we let U* := (/5*,0) with < p* < pQ, then it is easy to check that 

Ai([/*) = -pV(p*)<A2(f/*) = 0, 
and the hypothesis (HI) is satisfied as: 

VX2m ■ n{u*) = -p'ip*) = -7(p*)-'(^ - -)-^-' < 0. 

P Po 



(4.24) 



(4.25) 



(4.26) 
(4.27) 



(4.28) 



(4.29) 



(4.30) 



(4.31) 



Theorem 1.1 can thus be applied to ( 4.23| ) near the equilibrium U* or U 



A Proof of Proposition |2.2 



Proposition belongs to the folklore of finite-dimensional control theory (see in partic- 
ular Fillipov [|l^]). Since we have not found the exact required formulation in the literature, 
we give the proof in details for the sake of completeness. 

We begin with a few notations. 
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Definition A.l. 'P^ab) ^ L°°{a,b]^^) is defined as the set consisting of all piecewise con- 
stant vector functions on {a,b). Next ^(^i,^ C '^(^b) is defined as the set consisting of all 
piecewise constant vector functions on {a,b) with at most one nontrivial component, i.e., 
/ = (/i, • • • , JnY^ £ '^{ab) ^^'^ there exist p G N, indices ii, • • • ,ip G {!,••• , N}, 

constants /j\ , • • • , ff^ G M and a = Iq < ti < ■ ■ ■ < tp = b such that 

fit) = fid,, yt€{ti^i,ti),l = l,--- ,p, (A.l) 

where ei, • • • , eat denote the standard basis of M.^ . 

Now we deduce the following statement. 

Proposition A.l. •>^(oi) dense in L°°(0, 1;]R^) mi/i respect to the weak-* topology, more 



precisely, for any f G L°°(0, 1;]R^), there exists a sequence {f^}'^i C ^^i^ such that 



lim 

k^oo 



[ f^{t)-h{t)dt= f f{t)-h{t)dt, V/i G l1(0,1;]R^), (A.2) 

II/'|Il-(o,1;M-) < C^II/IIl-(0,1;M-)' VA; G N. (A.3) 
Proof: It is classical that "P/q D is dense in L°°(0, 1;M^) for the weak-* topology (moreover 



one can require (A.3) to hold on an approximating sequence). Hence it suffices to prove that 
^[^i^ is dense in 'P^qi) with respect to weak-* topology. To do this, we first prove ( |A.2| ) in 
the simpler case where / is a constant function: 

/(t)=7=(7i,--- Jjv)*'eM^, VtG[0,l]. (A.4) 

For any A; G N, we let G -^(o i) be defined as 

7\t) := iV/.e., Vt . "-^;f + ' ).V; . {1, ■ ■ ■ , t}, i = 1, ■ ■ ■ A'. (A.5) 

Clearly converges weakly-* to / in L°°(0, 1; M^) as A; tends to oo. Now we treat the general 
case where / G V^^ ^-^ . We introduce times Q = t^ < ti < • • • < tp = 1 such that 



f{t) = f\ ytG{ti.i,ti),l = l,--- ,p, (A.6) 



where the / are constants. 

From the previous arguments, we can obtain by translation and scaling that there exists 
{/ }T=i -^{ti 1 ti) such that / converges weakly-* to /' in L°°(t;_i, t;; M^) as k tends to 
oo. Finally, for any A; G N, we let 

/'=(*):= /(t), tG{ti_,,ti),l = l,--- ,p. (A7) 

It is obvious that {f^}'^^i C ^^i) and f^ converges weakly-* to / in L°°{0, 1;M^) as k tends 
to oo, i.e., ([A^) holds. Observe that (pj) holds. □ 
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Back to the proof of Proposition 2^. Since a G L°°(0, 1;]R" the solution z G 
C'^([0, 1];IR") to the ordinary differential equation ( |2.43| ) is Lipschitz continuous, since 



z{s) 



Y,amrM^))dO, VsG[0,l]. 



(A. 



Let a be such that (|2.42| ) holds. By Proposition AA, there exists a sequence {a'^}^]^ C 
with the notation := (a^^,--- ,a^_i,a^_^_i, - ■ ■ ,a^), which converges weakly-* to a in 
L°°(0,l;M"-i) and 



ll" ||l°°(o,1;K"-^) - C'||a||Loo(o,i;Mn-i) < Ce, \/k G N. 
Let G C°([0, 5fc];M") be the solution to the Cauchy problem 



— = Y^a^is)rj{z'), zH0) = 0, 



(A.9) 



(A.IO) 



where G (0, 1]. By (A.9), z^ is uniformly Lipschitz continuous. 
If e is small enough, then by ( 2.42| ), we can deduce that 



S, = 1, 



(A.ll) 



that is, Zfc is defined on the whole time interval [0, 1], for all A; G N, and 

||/||vi/i.-(o,i;R") < C'e, V/cgN. (A.12) 

By the Arzela-Ascoli Theorem, there exists a subsequence {z'''^}'^^ C {^'^j^^ and z°° G 
C°([0, 1];M") such that z^' converges to z°° in C°([0, 1];M") as I tends to oo. Now it is 
straightforward to pass to the limit in ( [A.8| ) (even, the limit is unique). The conclusion 
follows. 
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